Vortex distribution in neutron stars: gravitational effects. 
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Abstract 

Neutron stars are supposed to be mainly formed by a neutron superfluid. 
The angular momentum is given by the vortex array within the fluid, and a 
good account of the observable effects is determined by its coupling with the 
crust. In this article we show that the gravitational field introduces important 
modifications in the vortex distribution and shape. The inertial frame drag- 
ging on the quantum fluid produces a decrease in the vortex density, which for 
realistic models is in the order of 15%. This effect is relevant for neutron star 
rotation models and can provide a good framework for checking the quantum 

effect of the frame dragging. 
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I. INTRODUCTION 



In most quantum systems the gravitational effects are extremely small, but there are some 
remarkable exceptions. One which has been already extensively studied is the generation 
of inhomogeneities in the inflationary phase of the universe, through the amplification and 
the transition to a classical regime of the quantum fluctuation of the inflaton field [jT|. In 
this case the relevance of the phenomena is related to the propagation of modes with a 
wavelength comparable to the radius of the horizon of the universe, i.e., the gravitational 
and the quantum characteristic lengths are of the same order. A natural question which 
arises is if there are other cases accessible from a phenomenological point of view, and 
the best candidate for an affirmative answer is provided by the neutron stars. There is a 
consensus not only about their existence but also about their direct relation with pulsars 0| . 
As it is well known, the pulsar dynamics shows a very interesting phenomenology closely 
related to neutron star models. Among the most suggestive phenomena displayed by pulsars 
are the glitches, which are closely linked to the internal dynamics of the star. The post-glitch 
relaxation of the angular velocity has a large time scale, which is interpreted as evidence of 
a crust containing a superfiuid medium Pj.The models usually assume a solid exterior crust, 
which is the directly observable zone, and inner regions mostly constituted by superfiuid 
neutrons. The superfiuid phase would be stable at the very high pressure in the interior of 
the star . This image is supported by theoretical computations as well as observational 
evidence Given that a superfiuid flow is irrotational, the angular momentum of the star 
must be supported by vortex lines. They are coupled to the exterior crust, and the different 
phenomena associated to the glitches are related with this coupling 0. These models give 
a qualitatively correct description for the glitches dynamics. In this article we show that 
gravitational effects could be a very relevant ingredient to this description. Although this 
article is concerned with neutron stars, our conclusions can be applied to the rotation of 
boson stars. 

In the following section we give a brief introduction to the subject of the rotation of a 
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superfluid star, using the weak field approximation for tlie gravitational field. This approxi- 
mation is useful for understanding the phenomena, but it is not suitable for studying effects 
that involve the gravitational field in the interior of the neutron stars because of their high 
mass density. For this reason in Section III we analyze the rotation of a superfluid star 
using a different and more adequate approach. We consider there the exact expression for 
the gravitational interaction, and expand the metric in powers of the angular velocity. This 
is a small parameter for a neutron star, and thus we can restrict the expansion to linear tems 
without a significant error. Finally, in Section IV we study in detail the vortex distribution 
in the star and the effects of the gravitational field on this distribution using a covariant 
description for the superfluid. 



II. WEAK FIELD APPROXIMATION TO THE SUPERFLUID STAR ROTATION 

To get an insight into the physical aspects of this problem we will discuss it in the first 
place using a weak field approximation. Thus we will assume a post-Newtonian metric p[: 

g/iiy = Vf^u + , h^y <^l, (1) 

— * 

where the gravitational field is described by a vectorial potential, hi = /loi, and a scalar one, 
= In this case the Hamiltonian for a non relativistic particle becomes Q: 

^J = (l + 3^)^P±^ + m^. (2) 

Im 

Except for the factor (1 -|- 3 0), which gives the red shift effect, there is a clear similarity 
between this Hamiltonian an the corresponding one for a particle in an electromagnetic field, 
where the electric charge e should be replaced by the mass m , the vector electromagnetic 
potential A by —h, and the electric potential by the Newtonian one. Extending the analogy 
one could suppose that a gravitational Meissner effect, the expulsion of the field V x /i, 
is present among the gravitational phenomena that take place in a superfluid star, as was 
proposed in Ref. ||TU|. However, we can argue that this is not the case on the basis of the 
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Mach principle. It tells us that the field source tries to impose its own rest frame to the 
other bodies through the gravitational interaction. This implies that the gravitational effect 
on the bodies has a paramagnetic nature independently of the substances they are made of, 
and thus there is no such Meissner effect. This point will be formalized later. 

Let us now analyze the slow rotating star in this approximation. The energy-momentum 
tensor of the fluid is T^'^ = {p + p)u^u'^ + pg^^ ■, where p is the pressure and p the energy 
density. Expanding this tensor up to flrst order in the velocity and using the weak fleld 
approximation to the Einstein equations in the harmonic gauge, the equations for the fleld 
h become 

l\h=lQixG[{p + p)v + ]^{?,p + p)h] , (3) 
V. /I = . (4) 

Based on the analogy with a superconductor, the Hamiltonian implies that V x (-1? + 
h) = 0. Furthermore, in a superfluid star we have a stationary rotating fluid where V.v = 0. 
These two conditions, together with Eq.@ give v = —h. Thus Eq.(^) leads to 

Ah= -SttG {p-p)h. (5) 

As p > p, the paramagnetic character mentioned above is clearly stated. A similar result 
was obtained in Ref. 0], but there is a difference because the authors do not take into 
account the term proportional to h on the right hand side of Eq.(0). In the present context 
this term has the same weight as the term proportional to v, since v = —h, and it should 
not be neglected. 

The study of the rotation of the star reduces to the analysis of Eq.(^, which formally 
can be considered as a Schrodinger equation for h with a potential well proportional to p—p. 
In absence of vortices the phases induced in the superfluid by the rotation and the fleld h 
cancel exactly, and a rotational state of the star is given by a bounded regular solution of 
this equation that nullifles h at inflnity. In general the existence of this solution implies a 
deep or long potential well that leads to an unstable star that will collapse to a black hole. 
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At least this result shows that the weak field approximation is not suitable for analyzing the 
possibility of a vortexless rotational state, because non linear effects are important. 



III. SUPERFLUID STAR ROTATION 

In this section we will analyze the superfluid star rotation on the basis of an expansion of 
the metric tensor g^t in powers of the angular velocity Q, which in the case of neutron stars 
can be considered as a small parameter. This last point allows us to restrict the expansion 



up to linear terms in . At first order on the angular velocity the metric takes the form 



ds' = -e"^dt' + e'^-dr' + r\de' + sin '9 dip') - 2r' sin '6 oodipdt, (6) 

where uj{r, 6) represents the angular velocity that takes a free falling object from infinity to 
the point (r, 9), and which corresponds to the local inertial frames rotation with respect to the 
fixed stars. The potentials $ and A are even functions of the angular velocity, and therefore 
in this approximation they are the non-perturbed functions of r that can be computed 
directly from the non-rotating stellar model. For stars at the end of the thermonuclear 
evolution they can be determined from the equation of state and the central pressure. The 
equations in this case are: 



m(r) = / Airpdr , (7) 
Jo 

dp (p + p) {m + Anr^p) 



dr r{r — 2m) 

d(j) m + Ai^r^p 



(8) 
(9) 



dr r{r — 2m) ' 

e^^ = . (10) 

r — 2m 

The remaining component of the metric, a;, is given by the Einstein equation corresponding 
to the component R^t of the Ricci tensor. To reduce this equation to a simple and suitable 
form, we can introduce the angular velocity of the system, which can be precisely charac- 
terized by = m'^/m* in terms of the fiuid velocity . This magnitude f2 represents the 
angular velocity measured by an observer at rest with respect to the fiuid. The minimal 



energy configuration, and in consequence the stable one, has fl constant, which restricts 
the problem to the case of a uniform rotation. If we assume that the constitutive matter 
behaves perfect fluid, T^t is: 

T^t = -r^sm^9{{p + p)Q - puj). (11) 

Besides, using the azimuthal symmetry of the rotating star, we can expand u as follows: 

^ir.e)^Y..'(r)(-—,^). (12) 

where Pi is the Legendre polynomial of degree I. With these ingredients the perturbed 
Einstein equation becomes 



u\rr +(- - A' - ) u;\r +' (- + $'- A' - ^iillle^^) 
r ^ ^ 2r 

= 167re2^(i(p + 3p) J - {p + p)n 5[) . (13) 

The only information which remains to be introduced refers to the superfluid state of 
the star matter. The superfluid is characterized by the curved space-time covariant gen- 
eralization of Vx ?)= 0, which is satisfied by a superfluid in absence of vortices in a flat 
space-time. According to this the quadrivelocity must satisfy 

e^^-^'e/ («,,fc - uu,) = , (14) 

where ^ is a time-like Killing vector. Due to the symmetries of the system this relation 
reduces to = 0, which we can rewrite as 

-^(u'^g^^ + u'gt^) = 0. (15) 

Substituting the metric components by the expressions which correspond to (^), we 
obtain that the quantity r'^{u'^ — u^uj) is independent of r and therefore null. Hence we have 

^ = 1] = a; . (16) 

We will reconsider this relation further on to include the presence of vortices. By introducing 
it now in Eq.(|T^ and using (|^ we obtain the differential equation satisfied by u in the 
superfluid star: 
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Uj\rr +(- - A' - ) Uj\r + ^ J = Q. (17) 

In this equation we can substitute uj = , and develop up to the hnear terms in the 
fields, thus recovering Eq.(|^). For A and $ regular at the origin, i.e., nonsingular stars, the 
only solution with regular geometry is = 0. Therefore there are no solutions for rotating 



superfiuid stars without vortices. A similar result was obtained in Ref. [0 when analyzing 
boson stars. In this case there is not an effective equation of state and in that paper the 
boson star structure was computed using the Klein-Gordon equation. Besides, the energy 
momentum tensor is not isotropic because the radial and tangential pressures are not equal 
in the general case. However, as in Eg . (|13|) , only the radial pressure is relevant and the 
resulting equation is again Eq . {\L7\} . Due to the scalar field coherence the equation i7 = is 
satisfied, and hence the energy momentum tensor is a function of a; only. 

IV. ROTATION IN PRESENCE OF VORTICES 

Because the superfiuid star rotation cannot be achieved as a perturbation of its funda- 
mental state, we are going to study the rotation in presence of vortices. In order to consider 
the vortex contribution to the superfiuid star dynamics in curved space, we will follow an 
approach analogous to the one proposed by Weinberg Jl4|. The neutron (and proton) field 
has a U(l) global symmetry: 

^(x) — >e'^^{x), (18) 

that leads to the baryonic number conservation. The superfiuidity phenomena is related 
to the formation of a condensate that spontaneously breaks this global symmetry. This is 
similar to the superconductivity effect, where there is a spontaneous breaking of the U(l) 
electromagnetic gauge symmetry to Z2. Based on the physical picture for the neutron star 
matter, where neutron pairs have nonvanishing expectation values, we will assume that the 
U(l) baryonic symmetry is broken to Z2, the subgroup of transformations with A = and 
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A = vr. The spontaneous breaking of the symmetry leads to the existence of a Nambu- 
Goldstone excitation with zero energy in the hmit of vanishing momentum. The group 
transformation acting on the Nambu-Goldstone boson is: 

— > (j){x) + A . (19) 

As parametrizes U(l)/Z2 , and (p + n are taken to be equivalent. The U(l) invariant 
density Lagrangian is a function of the derivatives of and the U(l) fixed neutron fields 
^, and the baryonic Noether current is given by = ^(f^ • The most general density 
Lagrangian allowed by the symmetries is a nonlocal function of the field, but the nonlocality 
extends over a range of the order of the penetration length of the superfluid. Given that we 
are interested in the macroscopic fluid motion, we will only consider the local terms in the 
density Lagrangian that effectively describe the long range behavior. Such terms must be 
scalars and should be constructed as a contraction of covariant quantities. The only possible 
factors are the gradient of , the metric (not the curvature tensor because of the equivalence 
principle) and a number of fixed tensors that characterize the field \1/ and must satisfy the 
requirements of spherical symmetry, since they are determined by the unperturbed star. 
Let us call them A, \^^\ Ajfj, ... As was argued in Ref. the existence of an equilibrium 



configuration with vanishing gradients rules out the linear terms in its derivatives . The 
quadratic terms must not vanish since the system has a spontaneous symmetry breaking. 
Therefore, the density Lagrangian can be expanded as: 

^ = lifW + 9{\) Ag) 0'>^ + ...)■ (20) 

We only need the quadratic terms because is zero for the static star, and therefore 
it is a first order quantity of the angular velocity for a rotating one. Hence we have = 
/(A) 0_^ + (7(A) A*^^)^ (jy^y. Due to the spherical symmetry the tensor A|fj has the same angular 
dependence as the metric. Thus A'^^-'^J is a function of r and we can write = n{r) ,^ = 
n{r) dip(j). The current component j^p can also be written for small velocities as nQ{r)u^, 
where no(r) = is the baryon number density in the fluid rest frame and 7 is the 

determinant of the spacial metric tensor From here we obtain 
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Tiiv) 2 

u<P = -^Tdf<P = ' (21) 

no[r) m*[r) 

where m* = rio/n is a scalar quantity of dimension one that depends only on the fluid 
conditions. This parameter can be interpreted as the effective mass of the quasiparticles 
from the equivalence principle, and for a neutron star it takes a value of order of the neutron 
mass [0] |]T^. The Nambu-Goldstone field should be compared with the double the 
phase of the Ginzburg-Landau wave function, which is the origin of the factor two in the 



preceding equation [|1J] 



On the other hand we have 



SJ = !^ = M!l±i^!»!!. (22) 



This expression becomes 



^ = 2 ■ 2n t + ^ = * 2 ■ 2n ^ + ^ • (23) 

sm t7 M* m* sm t/ 

when we keep up to the first order terms in lo. The factor e* gives the redshift of the 
effective mass. From the rotation symmetry, the periodicity and the continuity of the field 
outside singularities we obtain = n(r, ^) where n(r, ff) is the sum of the indexes of the 
field singularities held by a circular closed path in a plane orthogonal to the rotation axis, 
centered on this axis, and containing the point (r, 6*). The n{r,6) function is also equal to 
the number of vortices surrounded by the closed path, because each vortex has a topological 
number one to minimize the energy. In consequence 

m* sm ^6 

This equation generalizes the expression (|16]) by including the presence of vortices. It also 
generalizes the formula describing the angular velocity and vortex density relation in fiat 
space. From this point of view the new ingredient is the second term on the left hand 
side of the equation, which corresponds to the dragging of the inertial frames due to the 
gravitational field. This relation shows that the phase introduced by the vortices in the 



superfluid is equal to the sum of the kinetic phase plus a gravitational phase. This phase is 
due to the metric that makes the covariant and contravariant components different, unlike 
the electromagnetic case and what can be interpreted from the weak field approximation 
where it is introduced by a connection. 

Given the angular velocity Q and the star structure at rest, we can state the vortex 
distribution in a superfluid star. First we compute u from Eq.(|l3D, with the reasonable sup- 
position that the vortex interaction energy is negligible with respect to the rotation energy. 
Once u is obtained we can compute the vortex distribution and shape from Eq.(^). The 
resulting vortex distribution minimizes the total energy of the fluid, leading to a macroscopic 
motion equal to the corresponding one without vortices. Since is a decreasing function 
of the radius, the decrease in the vortex density with respect to a fiat space is greater at 
the axis of the star. Besides this, the gravitational field produces a line vortex diffraction in 
such a way that these lines are not parallel to the rotation axis, as in a fiat space, except on 
the equatorial plane. The angle (3 between the rotation axis of the star and the direction of 
the vortex lines at the point (r, 6) is given by 

(l + /t(r)(/t(r) -2) sin2(0)J 
where K(r) = ^ (r^^ - + 

^ ' 2 \il—u: dr m* dr dr J 

Outside the star we have uj{r) = , where J is the total angular momentum. From 
here and applying Eq.(^) with = ^ and r = R, we get the total number of vortices 

N = m*{R) fnR^-'^) e-*(^). (26) 



Even if we do not know the star mass distribution, which is model dependent, it is possible 
to give an upper bound for the average number of vortices Ug per invariant area unit, dA = 
2iTre^dr. Given that e^^^-* > e"*^^-* for r < R. this upper bound is: 

Thus the relative decrease of the average vortex number density with respect to the one at 
a flat space, = ^^-^ ^ is: 
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where / is the star moment of inertia. To clarify the meaning of the formula we can put 
I = Mf^, with f the radius of gyration, and thus we get the decrease rate bound as 
, where Rg = ^-^^ is the Schwarzchild radius of the star. For different neutron star models 
quoted in Ref.( jT^)? this lower bound varies from 11% to 15%. The vortex density decrease 
is greater as the star is more relativistic. These results are interesting at least for two 
reasons. On the one hand, they can be a significant ingredient for analyzing the dynamics 
of neutron stars, provided the development of the high pressure nuclear matter theory and 
pulsar models are accurate enough for describing observations. On the other hand, they 
provide a framework for the verification of the equivalence principle at the level of quantum 
systems by checking the effect of frame dragging. 

V. ACKNOWLEDGMENTS 

We are grateful to G. Zemba for very interesting and illuminating discussions. This work 
was realized with a partial support from the Consejo Nacional de Investigaciones Cientificas 
y Tecnicas (CONICET), Argentina. 



11 



REFERENCES 



[1] E.W. Kolb and M.S. Turner, The Early Universe (Addison- Wesley Publising Co., Read- 
ing, Mass., 1990) 

[2] S.L. Shapiro and S.A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars: The 
Physics of Compact Objects (J. Wiley and Sons, New York, 1983). 

[3] M.A. Alpar, S.A. Langer and J.A. Sauls, Ap. J. 282, 533 (1984). 

[4] D. Pines and M.A. Alpar, Inside Neutrons Stars, a 1990 Perspective, in The Structure 
and Evolution of Neutron Stars, ed. D. Pines, R. Tamagaki, S. Tsutura (Addison- Wesley 
Publising Co., Reading, Mass., 1993). 

[5] M.A. Alpar, Superfiuid dynamics in neutron stars, in Neutron Stars: theory and obser- 
vation, ed. J. Ventura y D. Pines, 49 (Kluwer Academic Publishers, The Netherlands, 
1991). 

[6] D. Pines, M.A. Alpar, Nature 316, 27 (1985); G. Baym, C. Pethick y D. Pines, Nature 
224, 673 (1969); B. Link y R. Epstein, Thermally-driven neutron star glitches, |astro- 
ph/950802T| 



[7] R.I. Esptein, B. Link and G. Baym, Superfiuid Dynamics in the Inner Crust of Neutron 
Stars, in The Structure and Evolution of Neutron Stars, ed. D. Pines, R. Tamagaki, S. 
Tsutura (Addison- Wesley Pubhsing Co., Reading, Mass., 1993). 

[8] S. Weinberg, Gravitation and Cosmology (John Wiley, New York, 1972). 

[9] H. Casini y R. Montemayor, Phys. Rev. D50, 7425 (1994). 
[10] H. Peng, Gen. Rel. Grav. 15, 725 (1983); 22, 609 (1990). 
[11] C. Ciubotariu y M. Agop, Gen. Rel. Grav. 23, 405 (1996). 
[12] J. Hartle, Ap. J. 150, 1005 (1967); J.B.Hartle and K.S.Thorne, 



12 



[13] Y. Kobayashi, M. Kasai y T. Futamase, Phys. Rev. D50, 7721 (1994). 
[14] S.Weinberg, Prog.Theor.Phys.Sup., 43( 1986). 

[15] L.D. Landau and E.Lifchitz, The Classical Theory of Fields (Pergamon Press, Oxford, 
1980). 

[16] J. Wambach, T.L. Ainsworth and D. Pines, in Neutron Stars: theory and observation, 
ed. J. Ventura y D. Pines, 37 (Kluwer Academic Publishers, The Netherlands, 1991). 

[17] R.B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C38 , 1010 (1988). 



13 



